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Recently, stellar collapse involving black hole formation from massive stars is suggested to emit 
enormous fluxes of neutrinos on par with ordinary core-collapse supernovae. We investigate their 
detectability for the currently operating neutrino detector, SuperKamiokande. Neutrino oscillation 
... is also taken into account for the evaluation. We find that the event number is larger than or 

OO ' comparable to that of supernova neutrinos and, hence, black hole formation is also a candidate for 

, neutrino astronomy. Moreover, we find that the event number depends dominantly on the equation 

' of state used in the computations of the black hole formation. This fact implies that the detection 

of neutrinos emitted from the black hole progenitors is very valuable to probe the properties of the 
■4-^ , equation of state for hot and/or dense matter. 
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PACS numbers: 97.60.-s, 14.60.Pq, 95.85.Ry, 26.50. 



r-| ; I. INTRODUCTION 

Q , Detections of neutrinos from SN1987A [J, [2| have declared that neutrino astronomy is practical to study high cn- 
^ • ergy astrophysical phenomena. Since future detections of supernova neutrinos will bring useful information not only 
, qualitatively but also quantitatively, there are many studies on the estimations of the event numbers. In particular, 
about 10,000 events will be detected by SuperKamiokande (SK), which is currently operating neutrino detector, if a 
supernova occurs near the Galactic center now Q . ft is noted that the neutrino oscillation is discovered [1] and con- 
firmed by the solar, atmospheric, reactor and accelerator neutrinos [sl]. Since supernova neutrinos propagate through 
^ ' the stellar envelope and the earth, where neutrino flavor conversion occurs by the Mikheyev-Smirnov-Wolfenstein 
: (MSW) effect d, Q , we should take into account it to estimate the event numbers d S EHi ITsI T here are 
' studies taking into account even the resonant spin- flavor conversion the neutrino self-interaction [l5j and the 
, effects of the shock propagation [l^ for supernova neutrinos. The effects of neutrino oscillation on supernova relic 

neutrino background is also evaluated (l7| . 
_ It has recently been realized that neutrinos are emitted not only from iron-core collapse supernovae but also from 
rocygen- neon-magnesium (ONeMg) supernovae [ll], type la supernovae [l^ and black hole forming failed supernovae 
I [20I I21I [2^ . While the flavor conversions of neutrinos from ONeMg supernovae [1^ [13] and type la supernovae [l^ 
have been already studied, such a study has not been done of black hole progenitors (hereafter, we call them "failed 
supernova neutrinos") has not yet done. 

Stars more massive than ~ 25 solar masses (Mq) will form a black hole at the end of their evolution. Observationally 
^ these black hole progenitors have two branches, namely a hypernova branch and a faint supernova branch ,25:] . This 
5^ \ trend may suggest us that the strongly rotating massive stars result in the hypernova branch while non-rotating 
massive stars do the faint supernova branch. Stars belonging to the hypernova branch have large kinetic energy 
and ejected ^^Ni mass and they are implied to associate with the gamma-ray bursts. The numerical studies on the 
collapse of such stars need fully general relativistic computations for non-spherical situation and they are done merely 
for the simplified models so far 1261]. On the other hand, the kinetic energy and ejected ^^Ni mass are small for stars 
belonging to the faint supernova branch. The fates of non-rotating massive stars are classified also by a numerical 
study [27] as follows, a) Stars with main sequence mass, M < 2bMQ make explosions and produce neutron stars, b) 
Stars ranging 25Mq ^ M < AOAIq also result in explosions but produce black holes via fallback, c) For M > AOMq, 
the shock produced at the bounce can neither propagate out of the core nor make explosions. Faint supernovae may 
correspond to class b) of this classification. More recently, a black hole candidate whose mass is estimated as 24-33M0 
is discovered |28|] and it may be a remnant of class c). 
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In this study, we focus on the failed supernova neutrinos from class c) objects. In Refs.[2ll[22l. collapses of stars with 
40 and SOMq are computed by fully general relativistic neutrino radiation hydrodynamics under spherical symmetry. 
Having adopted the results of evolutionary calculations for the progenitors [H, [sO, [HI, [s^l as their initial models, 
they have evaluated the time evolution of emitted neutrino spectra as well as dynamics. As a result, they have 
concluded that the neutrino luminosity reaches ^ 10^^ erg/s, which is larger than those of ordinary supernovae [33| . 
Meanwhile these objects may be observed optically as the disappearance of supergiants because they do not make 
supernova explosions. In fact, a survey monitoring ^ 10^ supergiants is proposed very recently (3^ . If one supergiant 
is confirmed to vanish from this survey, corresponding data of neutrino detectors such like SK should be checked 
whether the burst of failed supernova neutrinos is detected. The signal of failed supernova neutrinos will provide us 
various valuable information, for instance, the equation of states (EOS) at high density and temperature As 
mentioned already, neutrinos undergo flavor conversion before the arrival at the detectors. One has to evaluate the 
oscillation effects of failed supernova neutrinos. 

We analyze the oscillation and observational aspects of failed supernova neutrinos using the results of Refs. [2l| 
(hereafter Paper I) and [IS (hereafter Paper II) . This paper becomes the first study to deal with them. It is well 
known that some parameters for the neutrino oscillation are not determined yet. Therefore, we also investigate their 
dependence. In order to deal the flavor conversion inside the stellar envelope, we utilize the results of evolutionary 
calculations for the progenitors, which are also adopted in Papers I and II. When neutrinos pass through the earth 
before detection, they undergo the flavor conversion also inside the earth. In this case, results of flavor conversions 
depend also on the nadir angle of the progenitor, which we examine its dependence. The event numbers of failed 
supernova neutrinos are evaluated for SK. 

We arrange this paper as follows. Issues on the neutrino oscillation are given in Sec. |TT1 In Sec. Illli we describe the 
methods to compute the event numbers at SK. The main results are shown in Sec. IIVI According to Papers I and II, 
the feature of neutrino emission differs depending on the EOS and progenitor model. Therefore the resultant neutrino 
number depends not only on the mixing parameters and nadir angle but also on the EOS and progenitor model. The 
main purpose of this study is to evaluate their dependence and ambiguities quantitatively. Sec. |V] is devoted to a 
summary. 



II. NEUTRINO OSCILLATION 



In this section, we introduce the general formulation of the flavor conversion by the MSW effect at first. Next, we 
evaluate the flavor conversion inside the stellar envelope using the results of evolutionary calculations for the stars 
with AOMq ^29., 30] and 50Mq [31|, [S^I- Here, a model from Ref. 29 is adopted as a reference model. Incidentally, 
this model is an initial condition of the numerical simulations in Paper I. Finally the earth effects are computed from 
realistic density profile of the earth [sll . 



A. General formulation of flavor conversion by MSW effect 

Here, we review the neutrino oscillation by the MSW effect in brief 0. Neutrino oscillation is caused by the 
discrepancy between the mass eigenstate and flavor eigenstate of neutrino. In general, the flavor eigenstate is related 
with the mass eigenstate as below, 
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where — sin^^j, — cos^y and Qij is the mixing angle for i,j = 1,2,3 (i < j). While sin^ ^12 « 0.32 and 
sin^ 023 ~ 0.5 are well measured from recent experiments, only the upper limit is given for sin^ ^13 < 2.0 x 10~^ 
0. Neutrino propagation is described by the Dirac equation for the mass eigenstate and we can take an extremely 
relativistic limit here owing to smaller masses of neutrinos comparing with their energy. Moreover, electron type 
neutrinos get the effective mass from the interaction with electrons when they propagate inside matter. Thus, the 
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time evolution equation of the neutrino wave functions can be written as, 
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where Gp, ndt) and E are the Fermi constant, the electron number density and the neutrino energy, respectively. For 
the anti- neutrino sector, the sign of rig (t) changes, t is an Affine parameter along the neutrino worldline. Am|j are the 
mass squared differences and they are experimentally determined as Am^i = 8 x 10~^ eV and |A?7i§]^| = 3 x 10"^ eV. 
Whether the sign of Am|]^ is plus (normal mass hierarchy) or minus (inverted mass hierarchy) is unclear under the 
current status [5|. Solving Eq. ^ along the neutrino propagation, we can follow the neutrino flavor eigenstate 

I 1 2 

and evaluate the survival probability of (/ = e, and their anti-particles) as p^f = \ {ipuf IV') | ■ 

A calculation of the survival probability can be divided to two parts, namely a probability that a neutrino generated 
as flavor / reaches to the earth with j-th mass eigenstate, P*{vf Vi), and a probability that a neutrino entering the 
earth with i-th mass eigenstate is detected as i^f, P®{vi Vf). Using them, the survival probabilities are expressed 
as. 



Y,P"{v,^v,)P®{v,^v^), (3a) 

3 

P*{^e ^ V^)P®{V^ ^ ^e), (3b) 



3 



i=l 
3 
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where represents the sum of and and Vx represents the sum of and Vt- In other words, py_^ is a probability 
that the neutrino generated as v^i is detected as i/^ or v^- on the earth. It is noted that we can assume ^-t symmetry 
owing to sin^ Oi^, « and sin^ 6*23 « 0.5. 



B. Flavor conversion inside stellar envelope 



Neutrinos emitted from the core of black hole progenitor propagate through the stellar envelope, which has the 
density gradient. In other words, the electron number density, ne(t), varies along the neutrino worldline. It is noted 
that the conversions occur mainly in the resonance layers. The matter density at the resonance layer is given as 

, / Am^\ /20McV\ /0.5\ , 9 ^ , ^ 

A-e.. ^ 6.6 X 10^ i^j-^j [~^^) [y^J ~ ^ ^ ' 

where Ye is the electron fraction. In the stellar envelope, there are two resonances, namely H resonance and L 
resonance, where H resonance lies at higher density (~ 2 x 10'^ g/cm'^) than L resonance does at lower density 
20 g/cm^). H resonance corresponds to the case of Am? = Am^^ and 6 = ^13 in Eq. (U), while L resonance does 
to An? = Am\i and 9 — 612. Thus, in case of the normal mass hierarchy, both resonances reside for the neutrino 
sector and the anti-neutrino sector has no resonance. On the other hand, in case of the inverted mass hierarchy, H 
resonance shifts to the anti-neutrino sector. If the density profile is gradual at the resonance layer, the resonance is 
"adiabatic" and the conversion occurs completely. When the resonance layer has a steep density profile, the resonance 
is "non-adiabatic" and the conversion does not occur. The flip probability of the mass eigenstates, Pp, is analytically 
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estimated as 



Pf = , (5a) 

— 1 



where r and rj-cs. are the radial coordinate and the position of the resonance layer, respectively. 

Using flip probabilities at H resonance {Ph) and L resonance (Pl), P*{vf Vi) in Eqs. ([3]) are calculated. Here 
we define matrixes P* and P* whose elements, P*j and P^j, are P*{vf Vi) and P*{vf i^i), respectively, and they 
are calculated as Q 

PlPh 1-Pl Pl{1-Ph) 
P* = I (f-PL)PH Pl (1-Pl)(1-Ph) |, (6a) 
1-Ph P„ 

0\ 

P* = I 1 , (6b) 
\0 ij 

in case of the normal mass hierarchy. On the other hand, in case of the inverted mass hierarchy, H resonance resides 
for the anti-neutrino sector and the survival probabilities are given as, 

/ Pl 1-Pl 0' 
P* = \ 1-Pl Pl I , (7a) 

(7b) 

In the absence of the earth effects, P®{vi vf) is displaced to |C^/i|^, where Ufi are elements of matrix U given 
in Eq. (fTbl). and Eq. ^ are rewritten as, 

P., « Ph [(l-sin2 0i2)PL +sin2 012(1 -Pl)] , (8a) 

p^- « l-sin2 6'i2, (8b) 

1 + Ph [(1 - sin^ 9i2)Pl + sin^ 612(1 - Pl)] 
P^. ~ 2 ' ^ 

_ 2 - sin^ 9,2 

P^-. ~ 2 ' ^ ' 

in case of the normal mass hierarchy and 

p,^ « (l-sin2 0l2)PL+sin2^^l2(l-PL), (9a) 
P.-^ « PH(l-sin2 0i2), (9b) 

V 
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1 + (1 - sin^ 0i2)Pl + sin^ ^12(1 - Pl) , 
« n ' (9c) 



in case of the inverted mass hierarchy. In Fig. [1] we show the results for the survival probabilities as a function of 
sin^ 6*13 for the reference model with E = 20 MeV. In this figure, the analytic formula, Eqs. ([8|) and ([9]), are compared 
with the numerical solutions of Eq. ^ by Runge-Kutta method. We can recognize that the analytic estimations well 
coincide with numerical results. This trend holds true for other progenitor models and neutrino energies, E, and we 
can safely use the analytic formula hereafter. In addition, L resonance is perfectly adiabatic (i.e., Pl ~ 0) for all of 
the progenitor models adopted here and for an energy range of our interest {E = 1-250 MeV). Thus, in case of the 
normal mass hierarchy p^^ « 0.68 and p^-^ « 0.84 are robust while p^^ « 0.32 and p^,^ « 0.66 are in the inverted mass 
hierarchy. 
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FIG. 1: Neutrino survival probabilities for the reference model with E — 20 MeV. Plots represent results of the numerical 
computations and lines do the analytic formula. 



C. Progenitor dependence 

In Paper II, the progenitor dependence of neutrino emission is investigated. Since the neutrino survival probabilities 
depend on the density profile of stellar envelope, we utilize the same stellar models of evolutionary calculations with 
Paper II. 40Mq model by Ref. [H (the reference model), 40Mq model by Ref. M and SOM© model by Refs. M, 113 
are named WW95, II95 and TUN07, respectively, in Paper II. We use the same notation in our paper. Incidentally, 
these models do not take into account the mass-loss whereas massive stars are suggested to lose their mass during 
the quasi-static evolutions. If strong a stellar wind sheds a large amount of the envelope, the density profile at the 
resonance layer will be affected. The mass-loss rate is thought to depend on the stellar metalicity, Z. Since TUN07 
is a model with Z = 0, the effects of mass-loss are absent. Incidentally the same model with TUN07 but Z = 0.02 
(solar metalicity) has also been computed In this paper, we adopt this model as TUN07Z in order to discuss its 
sensitivity. In addition, WW95 and H95 are solar metalicity models but without the mass loss. 

In Fig. [21 we show p,y^ and p^^ for the normal mass hierarchy and p,j^ and p^-^ for the inverted mass hierarchy as 
functions of sin^ 6*13 evaluated from the analytic formula. The neutrino energy is fixed to i? = 20 MeV for all models. 
We can see that the results of model WW95 differs from others. This is because the H resonance lies at the boundary 
of oxygen layer and carbon layer for model WW95 and its density gradient is steeper than those of other models. 

As already seen in Eq. the resonance layer depends on the neutrino energy, E. In Fig. [31 the energy dependences 
of p^^ for each progenitor model are shown for various values of sin^ 613 in case of the normal mass hierarchy, is 
insensitive to the progenitor model for sin^ 6*13 > 10""^ or sin^ ^13 < 10~^ whereas it is sensitive for 10~^ < sin^ 613 < 
10^^. In case of intermediate value of sin^ 613, there are energies which give a peak of Pi,^ for some models. Roughly 
speaking, H resonance of neutrinos with this energy lies at the boundary of layers. For instance, model TUN07 has 
the boundary of convective region and non-convective region near He shell burning layer. On the other hand, the 
envelope of model TUN07Z is already stripped and there are no such a boundary near the resonance layer. Thus, its 
result does not have a peak of p^^ even for the case of intermediate value of sin^ ^13 . These trends hold true for 
in case of the normal mass hierarchy, and for pir^ and Pi/- in case of the inverted mass hierarchy. 



D. Earth effects: nadir angular dependence 

In this paper, we evaluate the earth effects, P®{yi — > Uf), solving Eq. ||2]) along the neutrino worldline numerically 
by Runge-Kutta method using the realistic density profile of the earth [3^. It is noted that the density variation 
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along the neutrino worldline depends the nadir angle of the black hole progenitor. In Fig. Ul we show some examples 
for the results of the earth effects. We find that the spectral shape is deformed to wave-like and its shape varies with 
the nadir angle. This is because the typical length of neutrino oscillation becomes comparable to the size of the earth 
and the probabilities P®(yi —^Vf) become sensitive to the neutrino energy. The deformation of the spectrum occurs 
for most of the cases of the parameter sets. However, when the value of sin^ 6*13 is larger (> 10"'^) and H resonance is 
perfectly adiabatic, the deformation disappears for the neutrino sector in case of the normal mass hierarchy and for 
the anti-neutrino sector in case of the inverted mass hierarchy. Incidentally, these features of deformation are very 
similar with the case of the ordinary supernovae [l^, and ONeMg supernovae [l^l . 



III. SETUPS FOR FAILED SUPERNOVA NEUTRINO DETECTION 



A. Neutrino flux and spectrum before oscillation 



In this study, results in Paper I and II are adopted as the number spectra of the failed supernova neutrinos before 
the oscillation, dN^^/dE. The numerical code of general rclativistic t^-radiation hydrodynamics adopted in these 
references solves the Boltzmann equation for neutrinos together with Lagrangian hydrodynamics under spherical 
symmetry. Four species of neutrinos, namely i^e, ^e, i^n and Dfj., are considered assuming that i^t and Vt are the 
same as and D^, respectively. It is noted that the luminosities and spectra of and 9^ at the emission (before 
the oscillation) are almost identical because they have the same reactions and the difference of coupling constants 
is minor. 14 grid points are used for energy distribution of neutrinos. We interpolate the energy spectra linearly, 
conserving the total neutrino number. We calculate the probabilities of flavor conversion and event numbers at the 
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FIG. 3: Energy dependence of the survival probabilities of for each progenitor model. Dotted, solid, long-dashed, short- 
dashed and dot-dashed lines represent the cases of sin^ 613 = 10~^, 10"*, 10~^, 10~^ and 10~*, respectively. The normal mass 
hierarchy is posited here. 



detector. Taking into account the oscillation, spectra of neutrino flux at the detector, dF^^/dE, can be expressed as: 
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where i? is a distance from a black hole progenitor to the earth. In this study, we set i? = 10 kpc, which is typical 
length of our Galaxy. 

In Papers 1 and 11, the core collapse of massive stars with 40 and bOM^ [2^ [s^, [3l|, [s^] have been computed using 
two sets of supernova EOS's. They have concluded that the EOS at high density and temperature can be constructed 
from the duration of failed supernova neutrinos. In these references, they have adopted the EOS by Lattimer & Swesty 
(1991) (LS-EOS) 38] and EOS by Shen et al. (1998) (Shen-EOS) [39, 40]. LS-EOS is obtained by the extension of the 
compressible liquid model with three choices of the incompressibility. The case of 180 MeV, which is a representative 
of soft EOS, has been chosen in Paper 1 and 11. Shen-EOS is based on the relativistic mean field theory and a rather 
stiff EOS with the incompressibility of 281 MeV. Progenitor models adopted in Paper I and IT are already stated in 
Sec. Ill CI and their calculated models are listed in Table |T1 We use the same names for the models given in Paper 
1 and II hereafter. Since Shen-EOS is stiffer than LS-EOS, the maximum mass of the neutron star is larger and it 
takes longer time to form a black hole. Therefore the total energy of emitted neutrinos is larger for Shen-EOS than 
LS-EOS. In Fig. [51 time-integrated spectra without the neutrino oscillation are shown for models W40S and W40L. 
The plots for are not shown because they are almost the same with those for v^- However, this feature is not the 
case after the oscillation is taken into account. 
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FIG. 4: Energy dependence of the neutrino survival probabilities with the earth effects for the reference model. Solid, short- 
dashed, long-dashed and dot-dashed lines represent the cases of the nadir angle 0nad ~ 0°, 30°, 60° and 90°, respectively. Here 
Snad ~ 90° is equivalent to the case without the earth effects. 



TABLE I: Summary of calculated models in Paper I and II. For the names of progenitors and EOS's, see in Sec. Ill CI and 
Sec. nil Al respectively. Mprog and tsn are the mass of progenitor and the time at the black hole formation measured from 
the core bounce, respectively. _E^°' is the total energy of emitted neutrinos for all species. The averaged neutrino energy of all 
species is defined as (-Bail) = El,°^/Nl°^ , where Af*°* is the total number of emitted neutrinos for all species. 



Model 


Progenitor 


A/prog (Mq) 


EOS 


tBH (s) 


Er (ergs) 


(Sail) (MeV) 


W40S 


WW95 


40 


Shen 


1.35 


5.15 X 10"'' 


23.6 


W40L 


WW95 


40 


LS 


0.57 


2.03 X lO'^^ 


20.2 


T50S 


TUN07 


50 


Shen 


1.51 


4.94 X 10^^ 


23.7 


T50L 


TUN07 


50 


LS 


0.51 


1.66 X lO'^^ 


19.7 


H40L 


H95 


40 


LS 


0.36 


1.31 X 10^^ 


18.6 



B. Neutrino events at SuperKamiokande 

SuperKamiokande is a water Cherenkov detector, which is a descendant of Kamiokande, located at the site of 
Kamioka mine in Japan. The neutrino reactions in the detector which we take into account are as follows: 



+ p > e'^ + n, (ll^-) 

i^e + e > Ue + e, (lib) 

+ e > Pe + e, (11c) 

i^x +e — > i^x + e, (lid) 

i^x +e — > + e, (lie) 

— > e + i^F, (llf) 

ve + ^^O — ^ e+-f i^N, (Ug) 



We adopt the cross sections for the reaction (|llap from Ref. |4l|, for the reactions (jllf|l and ( |llgp from Ref. and 
for others from Ref. It is noted that the reaction (jllap makes dominant contribution to the total event number. 
Incidentally, the reaction (jllfj) has the largest cross section for neutrinos with the energy > 80 MeV although the 
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FIG. 5: Time-integrated spectra before the neutrino oscillation for models W40S (left) and W40L (right). Solid, dashed and 
dot-dashed lines represent the spectra of Ue., Ve. and Vx, respectively. 



ambiguities exist the theoretical cross sections of the neutrino-nucleus reactions. We assume that the fiducial volume 
is 22.5 kton and the trigger efficiency is 100% at 4.5 MeV and 0% at 2.9 MeV, which are the values at the end 
of SuperKamiokande \ The energy resolution was 14.2% for E^. = 10 MeV at this phase [ij and roughly 

proportional to yfEl ^10], where Eg, is the observed kinetic energy of scattered electrons and positrons. In this study, 
we choose an energy bin width of 1 MeV. 

IV. RESULTS AND DISCUSSIONS 

In this section, we report the time-integrated event number of failed supernova neutrinos at SK in the setup discussed 
above. The flavor conversions of neutrinos are computed from the density profiles of corresponding progenitors in 
Sec. Ill CI for 5 models listed in Table |T1 We investigate 2 additional models, Tz50S and Tz50L, whose models of 
emitted neutrinos are the same as T50S and T50L, respectively, but neutrino oscillations are calculated under the 
density profile of model TUN07Z. We note again that the progenitor model TUN07Z is the same as TUN07 but with 
the effects of the mass-loss. In Sec. IIV Al dependence on the undetermined parameters, namely sin 6*13 and the mass 
hierarchy, and the nadir angle ^nad is investigated. Dependence on the EOS's and progenitor models is shown in 
Sec. IIVBI Quantitative discussions on how these ambiguities affect the results are will be also made. 

A. Parameter dependence 

Here we focus on model W40S only since the qualitative features are similar among other models. In Fig.[6l we show 
the time-integrated event number of neutrinos for several sin^ ^13 values as a function of the nadir angle. The event 
number does not depend on the nadir angle very much and we will revisit later for the reason. It is also noted that 
the total event number is larger than that of the ordinary supernova (~ 10,000) Q because the total and averaged 
energies of emitted neutrinos are higher for the black hole progenitors. For instance, for the ordinary supernovae, the 
total and averaged energies of the neutrinos emitted as i7e are 4.7 x 10^^ ergs and 15.3 MeV, respectively [i^, while 
they are 1.42 x 10^^ ergs and 24.4 MeV in our model. Even the sum of total energies of emitted neutrinos over all 
species is lower in the ordinary supernova as 2.9 x 10^^ ergs. 

In case of the inverted mass hierarchy, the event number gets smaller for the larger sin^ ^13 . This is because the 
survival probability of Pe, which contributes to the event number most by the reaction (|llap . is almost zero for 
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FIG. 6: Nadir angular dependence of the time-integrated total event number of failed supernova neutrinos for the model W40S. 
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sin^ 013 > 10~^ (Fig.[2|). It is important to note that this tendency is different from the case of ordinary supernovae. 
For the ordinary supernova neutrinos, the flux of before the oscillation is larger than that of 9^ in the energy regime 
E > 25 MeV and the neutrino flux is not so small for this regime (see Fig. 35 in Ref. 3). For the case of inverted 
mass hierarchy and larger sin^ ^13, as seen in Fig. [2l half of Px converts to De and the event number of neutrinos with 
E > 25 MeV increases. Thus the resultant total event number of ordinary supernova neutrinos becomes larger. On 
the other hand, for failed supernova neutrinos, the neutrino flux is very small for the energy regime {E > 80 MeV) 
where the flux of Dx before the oscillation becomes larger than that of Pe (see Fig. [5]). 

The event number gets larger for the larger sin^ ^13 in the normal mass hierarchy while the difference is small. The 
difference owing to sin^ ^13 is mainly induced by the reaction (jllfjl . From Fig. [5l we can see that the differential 
number of ly^ before the oscillation is larger than that of i^e for E > 80 MeV. Thus the differential number of 
becomes larger by the oscillation for £^ > 80 MeV while it does smaller for < 80 MeV. Incidentally, the threshold 
of this reaction is high ~ 15.4 MeV and the cross section is roughly proportional to E^. Therefore the difference 
arises crucially for the high energy regime and resultant total event number is slightly larger for the larger sin^0i3. 
It should be remarked, however, that the precise evaluation for the high energy tail of neutrino spectra is difficult 
in the numerical computations and our results are thought to have some quantitative ambiguities. In addition, this 
feature is consistent with the case of ordinary supernovae. It is also noted that the mass hierarchy does not make any 
difference for sin^ ^13 < 10~^, as is the case of ordinary supernovae. 

In Fig. [3 we show the energy spectra of the time-integrated event number for several models. Results of the models 
with sin^ ^13 — 10~^ for the normal mass hierarchy and sin^ ^13 — 10~^ for the inverted mass hierarchy are not shown 
because they are very close to the result of the model with sin'^ 613 — 10~* for the normal mass hierarchy. The most 
outstanding case is the inverted mass hierarchy with sin^ ^13 = 10~^. The event number is much smaller than those 
of the other cases especially for the energy regime, -E < 50 MeV, and accordingly the spectrum becomes harder. The 
reason is again that the survival probability of Ve is almost zero and the original flux of is larger than that of i>x 
especially for this energy regime. This spectral feature may be useful for the restrictions of the mixing parameters, 
as already studied for ordinary supernovae 0]. 

In the case without the earth effects, spectra are not fluctuating while the cusps are artifacts due to the binning 
of original spectral data in Paper I and II. One exception is the model with sin^ ^13 = 10~^ for the inverted mass 
hierarchy, whose spectrum is fluctuating from 20 MeV to 30 MeV. This is because, as already mentioned in Sec. Ill CI 
the H resonance lies at the boundary of shells for the progenitor model WW95. The fluctuation originates from 
spikes in the upper left panel of Fig. [3l In the case with the earth effects for the nadir angle ^nad = 0°, 30° or 60°, 
fluctuating spectra are from the earth effects. This is because the spectral shapes of survival probabilities for each 
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FIG. 7: Energy spectra for the time-integrated event number of failed supernova neutrinos for the model W40S in the cases 
without the earth effects (upper left) and the nadir angle ^nad ~ 0° (upper right), 30° (lower left) and 60° (lower right). Solid 
and long-dashed lines represent the cases of sin'^ On = 10~* and 10~^, respectively, for the normal mass hierarchy. Dot-dashed 
and short-dashed lines represent the cases of sin^ di3 = 10~^ and 10~^, respectively, for the inverted mass hierarchy. 



neutrino species are deformed to wave-like as shown in Sec. Ill Dl In the case of the inverted mass hierarchy with 
sin^ 013 — 10~^, where there is no earth effect on the anti-neutrino sector, the deformation of spectra is not clear 
comparing with the other cases. Integrating over the neutrino energy, the fluctuation is smoothed out. Therefore the 
total event number does not depend on the nadir angle very much. 
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FIG. 8: Time-integrated total event number of failed supernova neutrinos for the normal mass hierarchy (left) and the inverted 
mass hierarchy (right). Error bars represents the upper and lower limits owing to the different nadir angles. The upper and 
lower sets represent models W40S and W40L, respectively. 



B. EOS and stellar model dependence 



As found in Papers I and II, the luminosity and duration time of failed supernova neutrinos depend on EOS and 
progenitor model. In Fig. \8\ the time-integrated event numbers of models W40S and W40L are shown for various 
parameter sets. The total event number of model W40S is at least 29,124 while that of model W40L is at most 
16,779. This large difference arise from not only its long duration time but also its high average neutrino energy for 
model W40S. For all parameter sets, event number of model W40S is at least 2.5 times larger than that of model W40L. 
Estimating from the studies of postbounce evolutions (< 1 sec) of the neutrino emission [ssj and protoneutron star 
coolings [46], the difference by the EOS is not so large for ordinary supernovae. Unfortunately, there is no previous 
investigation of an EOS dependence of supernova neutrinos from numerical computations of the successful explosion 
and the successive long-term (> 10 sec after the bounce) emission. However, it is likely that the black hole formation 
is a better site to probe the EOS of nuclear matter. 

In Table HIl the summary of event numbers for all models are shown. We can see that the event number of failed 
supernova neutrinos is mainly determined by properties of the EOS. Especially for the models by Shen-EOS (W40S 
and T50S), the difference is very minor. As discussed in Paper II, the duration time of the neutrino emission, which 
is mainly determined by the EOS, is also affected by the mass accretion rate after the bounce. The accretion rate 
depends on the density profile of a progenitor and T50S is reported to have lower accretion rate than W40S for late 
phase. If the accretion rate is lower, the neutrino duration time becomes longer but the neutrino luminosity is lower. 
This is because the neutrino luminosity is roughly proportional to the accretion luminosity. Therefore both effects 
are canceled out and the resultant event number becomes similar for models W40S and T50S. In case of LS-EOS 
(W40L, T50L and H40L), the event numbers differ by a factor of 1.8 at most and they are roughly proportional to 
the neutrino duration time. For these cases, the black hole formation occurs fast (< 0.6 s) after the bounce and the 
difference in neutrino luminosities induced by the accretion rate is not so large for this phase. Differences between 
the progenitor models are smaller than those between the EOS's for the models investigated in this study. 

We show in Fig. [5] the spectra of each model for several parameter sets. We can see the clear difference between two 
EOS's. Moreover, In the case of the inverted mass hierarchy with sin^ ^13 = 10~^, another interesting feature is seen 
in the spectra (lower left panel of Fig. [9|). It is notable that peaks are mismatched in the comparison between the 
spectra of models W40S and T50S. This comes from the different dependence of the neutrino survival probabilities 
inside the stellar envelope as shown in Fig. [2l The density at H resonance for neutrinos with the energy which gives the 
peak of spectrum in Fig. [S] roughly corresponds to that of the boundary of layers, where the density profile becomes 
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TABLE II: Time-integrated event number of failed supernova neutrinos for all models with various cases of the parameter sets. 
Two values of each column represent the upper and lower limits owing to the different nadir angles. 
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steep. 

In order to study the impact of the mass-loss of a progenitor, we compare models T50S, Tz50S, T50L and Tz50L 
in Fig. [ini Results for models T50S and Tz50S are the same except the case of the inverted mass hierarchy with 
sin^ ^13 = 10^'^. This is also the case for the comparison of models T50L and Tz50L. For the regime of a neutrino 
energy E = 30-40 MeV, the spectrum of model T50S has a bump, which does not appear for model Tz50S. This 
bump comes again from a steep density profile due to the onion-skin-like structure of the progenitor. The envelope of 
progenitor model TUN07Z is stripped by the effects of mass-loss and corresponding boundary is absent. Thus, only 
when the mass hierarchy is inverted and sin^ ^13 ~ 10~^-10~'', we have chances of probing the structure of a stellar 
envelope and effects of mass-loss. Incidentally, mass-loss may affect also the structure of more inner region. If this is 
the case, the dynamics and neutrino emission themselves may be changed. 

In conclusion, the event number of failed supernova neutrinos depends primarily on EOS. Ambiguities on the mixing 
parameters also affect the result and the spectra of detected neutrinos may be useful for the restrictions of them, 
while they can also be determined from other experiments. The event number reflects the properties of progenitor 
models though its dependence is weaker than that of EOS. For the limited case of the mixing parameters, we may 
be able to probe the structure of a stellar envelope. There is weak dependence on the nadir angle, which should be 
determined from the identification of the progenitor astronomically. Incidentally, Papers I and II discussed mainly 
about the time evolutions and duration times of the neutrino emission. Here, we have focused on the features of 
time-integrated event numbers and spectra, which are complementary to the previous studies. 

We remark that the conclusion given above is based on the results of some restrictive models and more detailed 
investigations are necessary to examine how general our conclusion is. Incidentally, EOS including the new degree of 
freedom such as hyperons |47'| or quarks [isj are shown to have different features at hot and dense regime recently. A 
new EOS of the nuclear matter for astrophysical simulations is also proposed to construct based on the many body 
theory, which is a different approach from Shen-EOS and LS-EOS [49]. However, we stress that the ambiguities and 
model dependences which can be taken into account under the current status have been fully involved and we feel 
that the essential point has been clarified in our investigations. 

Finally, we speculate the observational feasibility. The event number of models computed under Shen-EOS is at 
least 29,004 while that of LS-EOS is at most 16,779. Let us assume that the signal of failed supernova neutrinos 
is discovered from the data of the detector and its progenitor is identified in future. These facts indicate that the 
difference of EOS can be distinguished even if the ambiguities on the mixing parameters and progenitor models are 
not constrained well. If the mixing parameters are determined till that time by another way, the discrimination of 
EOS is more promising. In this analysis, the determination of the distance from the progenitor is a key factor. If we 
are lucky, the progenitor will have already been monitored by the surveys as presented in Ref. [s^ . Even if there are 
no ex-ante surveys, we can determine the direction of the progenitor to some extent by the neutrino detection itself 
[50| . Then the progenitor may be identified from old photographs. In addition, the nadir angle of the progenitor can 
also be determined so that we can remove the ambiguities of the nadir angle. 

The rate of these black hole forming collapses is highly uncertain but will be less than that of ordinary supernovae. 
Assuming the Salpeter's initial mass function, the black hole formation rate is ~ 25% that of ordinary supernovae [l^l ■ 
Since the supernova rate of our galaxy is ~ 0.03 /yr [5l|, we have to wait ~ 130 years for the black hole formation. 
One solution to this rate problem is to observe the nearby galaxies by the future detectors with larger fiducial volume. 
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FIG. 9: Time-integrated spectra for the total event number of failed supernova neutrinos for the normal mass hierarchy with 
sin^ 013 = 10~* (upper left), the normal mass hierarchy with sin^ S13 — lO"'^ (upper right), the inverted mass hierarchy with 
sin^ 6^13 = 10~^ (lower left) and the inverted mass hierarchy with sin'^ ^13 = 10~^ (lower right). Results obtained without the 
earth effects are shown. Solid, long-dashed, short-dashed, dot-dashed and dotted lines represent models W40S, W40L, T50S, 
T50L and H40L, respectively. 



Andromeda galaxy (M31) lies 780 kpc away. In the most optimistic case of Shen-EOS, the event number of neutrinos 
emitted from it is roughly estimated as ~ 8.2 for SK. If we consider a future detector with a scale of 5 Mton (e.g. 
the proposed Deep-TITAND detector jH, HI]), the event number becomes ~ 1900. This will be enough number to 
compare the LS-EOS case with ~ 500 events. Since the cumulative supernova rate within the distance to M31 is 
~ 0.08 /yr [5l|, we should wait ~ 50 years. If we extend further out to 4 Mpc, the event number becomes 70 for the 
most optimistic case. Since the supernova rate within 4 Mpc is 0.3 /yr, we should wait 13 years for the black 
hole formation. 
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FIG. 10: Same with Fig.[7]but for models T50S (upper left), Tz50S (upper right), T50L (lower left) and Tz50L (lower right) 
without the earth effects. 



V. SUMMARY 



In this study, we have evaluated the event number of neutrinos emitted from the stellar collapse involving black 
hole formation (failed supernova neutrinos) for the currently operating neutrino detector, SuperKamiokande. Since 
features of emitted neutrinos are affected by the property of equation of state (EOS) and the progenitor model f2ll.[2^. 
we have investigated the impacts of them under the currently available models. We stress that this is the first study to 
take into account the effects of flavor mixing for failed supernova neutrinos. Neutrino oscillations have been computed 
using not only the realistic envelope models of progenitors but also the density profile of the earth. We have also 
examined the dependence on the mass hierarchy and the mixing angle, sin^ ^13, in the range allowed up to now. 

As a result, we have found that the total event numbers of failed supernova neutrinos are larger than or comparable 
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to those of supernova neutrinos for all cases computed in this study. This result indicates that the black hole 
formation induced by the stellar collapse is the candidate for the neutrino astronomy as well as the ordinary core- 
collapse supernovas in spite of the shorter duration time of the neutrino emission. It is also important to note that 
the failed supernova neutrinos may hardly be negligible to estimate the flux of relic neutrino backgrounds. We have 
also found that the event number of failed supernova neutrinos depends primarily on EOS. Ambiguities on the mixing 
parameters also affect the result especially for the case of the inverted mass hierarchy and larger sin^ 6*13. Modifications 
by the progenitor models are not so drastic comparing with those by EOS, however, the structure of their envelopes 
affects the result for the limited case, the inverted mass hierarchy and sin^ ^13 ~ 10~^-10~^. While spectra of events 
are deformed by the earth effects, their impacts arc minor for the total event numbers. These results imply that failed 
supernova neutrinos are favorable to probe the properties of EOS for hot and/or dense matter. 
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